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Abstract. The D-dimensional Smorodinsky- Winternitz system, proposed some years ago 
^ ' by Evans, is re-examined from an algebraic viewpoint. It is shown to possess a potential 

algebra, as well as a dynamical potential one, in addition to its known symmetry and 
dynamical algebras. The first two are obtained in hyperspherical coordinates by introdu- 
, cing D auxiliary continuous variables and by reducing a 2D-dimensional harmonic oscillator 

Hamiltonian. The su(2Z?) symmetry and w{2D) ©s sp(4£',]R) dynamical algebras of this 
'(—( ', Hamiltonian are then transformed into the searched for potential and dynamical potential 

^ i' algebras of the Smorodinsky- Winternitz system. The action of generators on wavcfunctions 

, is given in explicit form for D ^ 2. 

' Key words: Schrodingcr equation; superintegr ability; potential algebras; dynamical potential 

algebras 

2010 Mathematics Subject Classification: 20C35; 81R05; 81R12 

> 

0\ ■ 1 Introduction 

(N 
O 



In classical mechanics, a Hamiltonian H with D degrees of freedom is said to be completely 
integrable if it allows D integrals of motion X^, i^i = 1,2, D, that are well-defined functions 
on phase space, are in involution and are functionally independent (see, e.g., [1]). These include 
the Hamiltonian, so that we may assume Xd = H. The system is superintegrable if there exist k 
additional integrals of motion Yy, v = 1,2, . . . ,k, 1 < k < D — 1, that are also well-defined 
functions on phase space and are such that the integrals H, Xi,X2, ■ ■ ■ , X^-i, Yi,Y2, . . . ,Yii. are 



^ , functionally independent. The cases k = 1 and k = D — 1 correspond to minimal and maximal 



superintegrability, respectively. 

Similar definitions apply in quantum mechanics with Poisson brackets replaced by commuta- 
tors, but H, X^, and Y^, must now be well-defined operators forming an algebraically indepen- 
dent set. Maximally superintegrable quantum systems appear in many domains of physics, such 
as condensed matter as well as atomic, molecular, and nuclear physics. They have a lot of nice 
properties: they can be exactly (or quasi-exactly) solved, they are often separable in several co- 
ordinate systems and their spectrum presents some "accidental" degeneracies, i.e., degeneracies 
that do not follow from the geometrical symmetries of the problem. 

The most familiar examples of such systems are the Kepler-Coulomb [21 [3lll] and the oscilla- 
tor [H [6] ones. Other well-known instances are those resulting from the first systematic search 
for superintegrable Hamiltonians on E2 carried out by Smorodinsky, Winternitz, and collabora- 
tors O [El [9] and from its continuation by Evans on E3 [lO] . These studies were restricted to those 
cases where the integrals of motion are first- or second-order polynomials in the momenta. Later 



*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S**)". The full collection is available at |http:/ /www.emis.de /journals/SlGM A/S4.html , 
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on, many efforts have been devoted to arriving at a complete classification of these so-called 
second-order superintegrable systems (see, e.g., [T H [T ^ [TBt [T^ \W \ [TB I fT7]). 

Only recently, the pioneering work of Drach |18[ [T9] on two-dimensional Hamiltonian systems 
with third-order integrals of motion has been continued |20[ [2T] . Nowadays the search for D- 
dimensional superintegrable systems with higher-order integrals of motion has become a very 
active field of research (see, e.g., [SgiSllHESlESlETllMlESlEQlEIlE^ 

In the present paper, we plan to re-examine from an algebraic viewpoint one of the clas- 
sical examples of Z)-dimensional superintegrable quantum systems, namely the Smorodinsky- 
Winternitz (SW) one [3 [HI [U [lOl [33l IM] , which may be defined in Cartesian coordinates as 



Here w, /ci, /c2, • • • , ko are some constants, which we assume to be real and positive. 

Several distinct algebraic methods may be used in connection with superintegrable systems. 
One of them is based on the fact that the integrals of motion generate a nonlinear algebra 
closing at some order |35[ [36l I37| . It has been shown, for instance, that for two-dimensional 
second-order superintegrable systems with nondegenerate potential and the corresponding three- 
dimensional conformally flat systems, one gets a quadratic algebra closing at order 6 [HI [121 [131 
[T4lll5j. Its finite-dimensional unitary representations can be determined |38j by using a deformed 
parafermion oscillator realization j39[l40|. thereby allowing a calculation of the energy spectrum. 
This procedure can be extended to higher-order integrals of motion and to the corresponding 
higher-degree nonlinear algebras [25[ [26] . 

Superintegrable systems may also be related [25l [26] to systems studied in supersymmetric 
quantum mechanics [41[ [32] or higher-order supersymmetric quantum mechanics [431 [^H [35| [^6] 
WI\ \^8\ [39l [50] I51| . hence can be described in terms of either linear or nonlinear superalgebras. 
As a consequence, supersymmetry provides a convenient tool for generating superintegrable 
quantum systems with higher-order integrals of motion j52[ [53] . 

The concept of exact or quasi-exact solvability |54[ [55] 156] . based on the existence of an 
infinite flag of functionally linear spaces preserved by the Hamiltonian or only that of one of 
these spaces, appears to be related to finite-dimensional representations of some Lie algebras 
of first-order differential operators, such as sl(2,]R), sl(3,M), etc. Although different from the 
concept of superintegrability, it can be related to the latter for some superintegrable systems 
(see, e.g., |27 [ 1491 [57] ). It is worth noting, however, that some alternative definitions of exact 
and quasi-exact solvability have been proposed for some specific superintegrable systems in 
connection with multiseparability of the corresponding Schrodinger equation |58[ 159]. 

The accidental degeneracies appearing in the bound-state spectrum of superintegrable quan- 
tum systems may be understood in terms of a symmetry algebra, which is such that for any 
energy level the wavefunctions corresponding to degenerate states span the carrier space of one 
of its unitary irreducible representations (unirreps) [601 161] . The generators of this symme- 
try algebra, commuting with the Hamiltonian, are integrals of motion, which may assume a 
rather complicated form in terms of some basic ones due to the fact that linear algebras are 
often preferrecJll (note, however, that nonlinear algebras may also be considered [62]). A familiar 
example of this phenomenon is provided by the so(4) symmetry algebra of the three-dimensional 
Kepler-Coulomb problem [21 [3] [1]. Another one corresponds to the su(3) symmetry algebra of 
the three-dimensional SW system [34] (or, in general, su(-D) for the Z)-dimensional one). 

In some cases, the symmetry algebra can be enlarged to a spectrum generating algebra (also 
called dynamical algebra) by including some ladder operators, which are not integrals of motion 

'^It is worth observing here that this may be seen as the obverse of the approach used in 1121 1131 1141 1151 
[35] [36] [37l [38] , where the generators are the basic integrals of motion but the algebra turns out to be nonlinear. 
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but act as raising or lowering operators on the bound-state wavefunctions in such a way that all 
of them carry a single unirrep of the algebra |63[ I64[ [55]. For the three-dimensional SW system, 
it has been shown [33] to be given by the semidirect sum Lie algebra w(3) (Bs sp(6,M), where 
w(3) denotes a Weyl algebra (or, in general, by w{D) (Bs sp(2L',M) in D dimensions). 

For one-dimensional systems, three other types of Lie algebraic approaches have been ex- 
tensively studied. All of them rely on an embedding of the system into a higher-dimensional 
space by introducing some auxiliary continuous variables and on the subsequent reduction of 
the extended system to the initial one, a procedure also used in discussing superintegrability 
(see, e.g., |24j ) . They work for hierarchies of Hamiltonians, whose members correspond to the 
same potential but different quantized strengths. The simplest ones are the potential algebras 
\6()\ [671 [68] , whose unirrep carrier spaces are spanned by wavefunctions with the same energy, 
but different potential strengths. Larger algebras, which also contain some generators connec- 
ting wavefunctions with different energies, are called dynamical potential algebras [BHlEQlElll. 




Finally, a third kind of algebras, termed satellite algebras |74| 175) . have the property that there 
is a conserved quantity different from the energy. 

Up to now, only the first one of these Lie algebraic approaches, namely that of potential 
algebras, has been applied to some D-dimensional superintegrable systems [76l [77l [78l [791 113 



The purpose of the present paper is threefold: first to apply this technique to the D- 
dimensional SW Hamiltonian (jl.ip . second to present for the same the first construction of 
a dynamical potential algebra in more than one dimension, and third to show very explicitly the 
action of both the potential and dynamical potential algebra generators on the wavefunctions 
in the two-dimensional case. 

The paper is organized as follows. In Section [2] the solutions, as well as the symmetry 
and dynamical algebras, of a 2D-dimensional harmonic oscillator are obtained in a suitable 
orthogonal coordinate system. In Section[3l they are transformed into the solutions, as well as the 
potential and dynamical potential algebras, of the D-dimensional SW system in hyper spherical 
coordinates. The D = 2 case is then dealt with in detail in Section [H Finally, Section [5] contains 
the conclusion. 

2 2D-dimensional harmonic oscillator 

Let us consider a harmonic oscillator Hamiltonian 



in a 2D-dimensional space, whose Cartesian coordinates are denoted by X^, fi = 1,2, .. . ,2D. 
For our purposes, it is convenient to consider it in a different orthogonal coordinate system, 
which we will now proceed to introduce. 

2.1 Harmonic oscillator in variables R, 0i, 62, • • • , Ai, A2, . . . , Xd 

On making the change of variables 

Xi = ii sin ^1 sin02 • • • sin0£i_i sin Ai, X2 = RsinOi sin02 ■ ■ ■ smO^^i cos Ai, 

^Some authors prefer to use the terminology of dynamical algebra of the hierarchy instead of dynamical 
potential algebra and to employ discrete variables, related to the quantum numbers characterizing the system, 
instead of continuous auxiliary variables. In this way, they get discrete-differential realizations of the algebras |72j . 
Other authors favour the use of nonlinear superalgebras [73] . 




[8ll[82l[83]. 



2D 
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X2u~i = R sin 9i sin 02 • • • smOu^u cos 9o~u+i sin Ajy, = 2, 3, . . . , — 1, 

X2,y = i?sin 01 sin 02 • • • sin0£)„,y cos 0£)_;^+i cos A^, = 2, 3, . . . , D — 1, 

X2Z)_i = iicos 6i sin Ad, ^20 = -Rcos 9i cos Ad, 



(2.1) 



where < < 00, < < f , = 1, 2, . . . , L» - 1, and < A^ < 27r, = 1, 2, . . . , L», H"^" can 
be rewritten as 



-dl - ^^^^dR - dl + [{2D - 3) cot 01 - tan0i]9ei 



R 



i/=2 



sin^ 01 sin^ 02 • • • sin^ 0,y_i l 



+ 



1 



+ [{2D -2u-l) cot 0^ - tan0^]90^ 

D-l 



+ 



sin 01 sin 02 • • • sin 6d-i 
1 



1 



u=2 



sin 01 sin 02 • • • sin 0£)-i/ cos^ 0£)-;/+i 



5L 



2 fl, 



COS^ Ol 

and is clearly separable. 

In the corresponding Schrodinger equation 

^osc^osc^^^ 6>, A) = 6>, A) 

with = 0102 • • • 0_D-i and A = A1A2 • • • Ad, we may therefore write 



D 



z = R^ 



(2.2) 



(2.3) 



where 



at*°-(i?,0, A) = -pl,_,+iM/°-(i?,0, A) 

and pi,P2, ■ ■ ■ ,Pd G ^- The normalization constant M"^^ in (j2.3p will be determined in such 
a way that 



j dV |^°^"(i?,0,A)p = 1, 



where 



2D 



dV =^1 dX^ = R^^^^dR 



n( 



sm ( 



\2D-2v-l 



COS 0,/(i0j, 



D 



WdK 

\u=l 



(2.4) 



(2.5) 



As shown in the appendix, the angular part of wavefunctions (j2.3p can be written as 

0nP^W= n n = nin2---nD-i, P = PiP2---PD, 

0K,fe.)(0^) = (cos0,)'^''-i(sin0,)^'^-5pi;:''-5''^+^-"-i)(_cos20,), 
where ni, 71-2, . . . , n^-i G N, 

du = [Pu[ + \i ly = 1,2, . . . ,D - 1, 



(2.6) 



(2.7) 
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bu = 2n^+i+ 2ni,+2H h 2nD-i+ \Pu+i\ + \Pu+2\ H h \pd\ + 5, 1^ = 1,2, . . . , D - 2, 

bo-i = \pd\ + i (2.8) 

— cos 29iy) denotes a Jacobi polynomial [HI], while the radial part can be 

expressed as 

C^£{z) = z^L(^^+^-'\z)e"'^^ (2.9) 

in terms of a Laguerre polynomial [M]. Here G N, while j is defined by 

j = m + 712 H h nz)_i + i(|pi| + IP2I H h IpdI) (2.10) 

and may take nonnegative integer or half- integer values. 
The corresponding energy eigenvalues are given by 

EZ) = 2{2nr + 2j + D). (2.11) 

We therefore recover the well-known spectrum of the 2D-dimensional harmonic oscillator 

E°^^ = 2{N + D), N = 2nr + 2j = 0,1,2, ... , 

whose levels, completely characterized by A^, have a degeneracy equal to 

Finally, the normalization constant in equation ()2.3p can be easily determined from some 
well-known properties of Laguerre and Jacobi polynomials j84j and is given by 



I \ 1/2 



™ \7rD{nr + 2j + D-l)\ 

ny + ay-}^)\{ny + hy + D -V -^)\ 



^ ( nj{2n^ + a^ + b^ + D-u-l){n^ + a^ + bu + D-iy-2y. \ 



X 



2.2 Harmonic oscillator symmetry and dynamical algebras 

As it is well known [5l [6], to each of the oscillator levels specified by N there corresponds 
a symmetric unirrep [N] of its su(2Z)) symmetry algebra. The generators of the latter 

EfMu = E^i, — ■^^t^,u ^ Epp, = 1,2, ... , 2D, 

with 



[Efj,u,Ep>y>] — 5u,^iE^yi — 5pyE^iy, EJj^j^ — E^p, 
are most easily constructed in terms of bosonic creation and annihilation operators 

al = ^{Xp-dx,), a^ = ^{Xp + dx,), fi = 1,2, . . . ,2D, (2.13) 

from 

Epu = = "i"^ + (2-^^) 

The harmonic oscillator Hamiltonian turns out to be proportional to the first-order Casimir 
operator Ci of u(2D), 

H''"' = 2Ci = 2^Epp = 2£. (2.15) 
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In the coordinates (|2.ip chosen to describe the oscillator, the so(2D) subalgebra of su(2Z)), 
generated by 

Lfj,,/ = — i {E^u — Eyfj) = —i{E^i, — Ey^), (2-16) 
such that 

is explicitly reduced. Its unirreps are characterized by 2j, which runs over A^, — 2, . . . , (or 1) 
for a given A^. The remaining generators of su(2D) may be taken as 

T^u = E^y + Eyfj^. (2.17) 

The operators 

E>l,u = (^li<^l, = a^Uu (2.18) 

act as raising and lowering operators relating among themselves wavefunctions corresponding 
to even or odd values of A^. Together with £^^,y, they generate an sp(4L',M) Lie algebra, whose 
(nonvanishing) commutation relations are given by 

[E>iiv,D\iyi] = &^L,iJ.'Ei,'i, + 6^yE^'i, + 5yyE^i^ + S^yE^i^. 

To connect the wavefunctions with an even A^ value to those with an odd one, we have to use 
the bosonic creation and annihilation operators (|2.13p . which generate a Weyl algebra w{2D), 
specified by 

The whole set of operators {E^y, Dji^, D^y^a^^^, a^, /} then provides us with the harmonic oscil- 
lator dynamical algebra, which is the semidirect sum Lie algebra w(2D) ©g sp(4D,R), as shown 
by the remaining (nonvanishing) commutation relations 

To apply the symmetry and dynamical algebra generators to the oscillator wavefunctions ()2.3p 
written in the variables i?, 0, A, we have to express the creation and annihilation operators ajt, 
in such variables. This implies combining the transformation ()2.ip with the corresponding change 
for the partial differential operators 

= sinA.a^'^'^) + cos A,a(^'2), 
9x2, =cosA,.a('''^) -sinA^a^^'^), = 1, 2, . . . , D, (2.19) 

where 

c)(i'^) = sin Q\ sin 02 • • • sin ^d-iS/j 
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1 ^"^ 

+ — CSC 6i CSC 02 • • • CSC 0p_i COS 9p sin ^p+i sin 9p^2 • • • sin Oo-idg^ , 
-I 

p=i 



d'^"'^'^ = sin 6^1 sin 6^2 

+ — CSC 01 CSC 02 



• • • sin ^z).;, cos Oo^v+idR 



■ ■ ■ CSC 0p_i COS 0p sin 0p+i sin 0p+2 • 



— CSC 6*1 CSC 02 • • • CSC Od-1 

R 



and 



^ sin0D-iy+i96iQ_„+i, 
5(^.1) = cos 0i5r--^ sin 0i5,,, 

5*-^'^^ = — CSC 01 CSC 02 • • • CSC 0£)_i9ai , 

S^'^'^^ = CSC 01 CSC 02 • • • CSC Od-v sec Bo-v+idx^ , f 

ri 

5(^'') = ^sec0ia,, 



sin 0D_,^ cos 0D-y+idep 
z^ = 2,3, 1, 



= 2,3,...,Z)-1, 



We shah carry out this transformation expHcitly for D = 2 in Section [H 

3 Reduction of the 2D-dimensional harmonic oscillator 
to the D-dimensional SW system 

To go from the 2L'-dimensional harmonic oscillator Hamiltonian H°^^ to some extended SW 
Hamiltonian let us first transform the original Cartesian coordinates X^, fi = 1,2, .. . ,2D, 
into some new ones x^, fi = 1,2,..., 2D, such that 



xi = rsmcpi sin<p2 • • • sin(p/)_i, 

Xi, = rsin(f>i sin</)2 • • • sin <f>D-u cos(f)D-u+i, 

xd = r cos 01 , 

XD+v = K, = 1,2,. . . ,D, 



u = 2,3,...,D-l, 



and 



R = V^r, 



— Yv 



, u = l,2,...,D-l, 



IT 



0<r<oo, 0<(pu<-, u = 1,2,... ,D -I, < A,. < 27r, u = l,2,...,D. 

Here r, (j)i,(j)2, ■ ■ ■ , (j^D-i are hyperspherical coordinates in the D-dimensional subspace (xi, X2, 
. . . ,xd)- The volume element in the transformed 2D-dimensional space is given by 



2D 



dv = JJ dxp^ = ^dr 
^l=l 



D-l 



n(sin<A.)^-''-'^i<A. 



D 



II dX, 



\u=l 



On making next the change of function 
^'(r,0,A) = 0^/^'^'''''{R,e,\), 

with 



D-l 



O = (wr)^ J] (sin f/.^)^-'' cos ( 



(3.1) 



(3.2) 



(3.3) 
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the harmonic oscihator wavefunctions ^°^'^{R,9, X), hving in a Hilbert space with measure dV 
given in (|2.5|) . are mapped onto some functions ^(r, 0, A), hving in a Hilbert space with mea- 
sure dv defined in (j3.ip . As a consequence of (j2.4p . we obtain 



j du|^(r,0,A)|2 = 1. 



By this unitary transformation, the harmonic oscillator Hamiltonian H°^^ is changed into 
H/UJ = Ol/2^osc0-l/2 ^34) 

and similarly for other operators acting in the harmonic oscillator Hilbert space. A straightfor- 
ward calculation leads to the result 

H = -dl - ^^^dr -\\dl^ + {D-2) cot cPid^, 



D-1 ^ 

+ Yl ■ 2, ■ 2 , [^i + {D-i^-l) cot (l)ud^^] 

sm (pi sm 02 • • • sm (pu~i 

sin^ 4>i sin^ <?!>2 • • • sin^ (t>D-i '^^ 4^ 

^ sin2,/)isin2 02---sin2 0D-^cos2</.B_^+i V 4y cos^i V 4 
-fwV^ (3.5) 

The eigenvalues of H are directly obtained from (j2.1ip as 

En^j = 2uj{2nr + 2j + D), Ur = 0,1,2, ... , j = 0,i,l,|,.... (3.6) 

The corresponding wavefunctions can be derived from (12. 3j) . (12. 6p . (j2.7p . (|2.9j) . (I2.12p . (|3.2p . 
and ()3.3p and read 

*n.np(r,0,A) = AA„,„p4-^)(z)$(r)(0) I^He^P-^+^^'^j , 

= (^)'^^ L(?/+^-i)(.)e-^, z = ..r^, 

D-l 

1^=1 

= n(cos0.)'^''(sin,^,)^''+5(^— i)pi:^-^'^''+^-'^-i)(-cos2<^,), 
j/=i 

AT^^np = ^^^^AC-p, (3.7) 

where n = nin2 • • • no-i, p = P1P2 • ■ ■ PD, nr,ni,n2, . . . ,nD^i G N, pi,P2, ■ ■ ■ ,Pd G while j, 
and a,y, are defined in (I2.10p and ()2.8p . respectively. 

In the subspace of functions 'I'n,.np('^i 0, •^) with fixed p, the Hamiltonian H, defined in (j3.5p . 
has the same action as the D-dimensional Hamiltonian 

//(k) = - ^^dr -^ldl+{D-2) cot <Pid^, 
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D-l ^ 

+ V ■ 2 , -2^ [^l. + {D-iy-l) cot (j)^d^ 

^-^ sm (pi sm 02 • • • sm q>v-\ 



v=2 



D-l 



v sm" (/>! sm" (pi - ■ ■ sm" sin^ c/ii sin^ 02 " " " sin^ ^d-v cos^ ^D-i^+i 

+ ^-R + u?r^ 

cos^ 01 

where we have defined \i = kik2 ■ ■ ■ kjj and 



A:. = V^D-.+i - i' z^ = l,2,...,D. (3.8) 

The latter Hamiltonian is but the SW one (jl.ip . expressed in hyperspherical coordinates r, 
01, 025 • • • ) 0D-1- We conclude that if is an extension of H^^\ resulting from the introduction 
of D auxiliary continuous variables X^, = xd+p, v = 1,2, ... ,D, and that, conversely, is 
obtained from B. by projecting it down into the D-dimensional subspace (xi,2;2, . . . 

As a by-product of this reduction process, we have determined the wavefunctions '^^}niT^ <p) 
of H^^^ in hyperspherical coordinates. Equation (|3.7p may indeed be rewritten as 

D 

1^=1 

'^'nX{r,ct>) =J\fi':lzii\z)^^^\ct^), M^'^l = (27r)^/2AA„,„p, (3.9) 

with k and p related as in (|3.8p . 

By a transformation similar to (j3.4p , the generators E^i, of the harmonic oscillator symmetry 
algebra su(2L') are changed into some operators acting on ^'„^np(?', 0, A). Since the latter may 
change rir and j separately (provided their sum rir + j = N/2 is preserved), this means in 
particular (see equation (12.10p ) that the p^'s (hence the k^'s) may change too. The transformed 
su(2D) algebra may therefore connect among themselves some wavefunctions of H belonging 
to the same energy eigenvalue ()3.6p . but associated with different reduced Hamiltonians H^^\ 
We conclude that it provides us with a potential algebra for the SW system. Similarly, the 
transformed w{2D) (Bs sp(4Z),M) algebra will be a dynamical potential algebra for the same. 

As a final point, it is worth observing that in the harmonic oscillator wavefunctions (12. 3p . the 
quantum numbers P;/, = 1, 2, . . . , Z), run over Z. On the other hand, in (jl.ip . the parameters kiy, 
V = 1,2, . . . , D, have been assumed real and positive. From equation (13. Sp . however, it is clear 
that pd-u+1 = would lead to an imaginary value of k^ and to unphysical wavefunctions (|3.9p . 
while |pz)_i,+i| and —\pD~u+i\ with \p£)-u+i\ > 1 would give rise to the same k^, hence to some 
replicas of physical wavefunctions (j3.9p . The correspondence between the harmonic oscillator 
wavefunctions and the extended SW ones is therefore not one-to-one. This lack of bijectiveness 
is a known aspect of potential algebraic approaches (see [69] where this phenomenon was first 
pointed out). 



4 The two-dimensional case 

4.1 Harmonic oscillator symmetry and dynamical algebras 

To deal in detail with the two-dimensional case, it is appropriate to rewrite the four-dimensional 
harmonic oscillator wavefunctions ^'^^,^p_^ ,p^{R,6,Xi,X2) (with j = n + ^{\pi\ + IP2I)) in an 

^Strictly speaking, this is true only for those k^'s that can be written in the form l|3.8|) with integer 
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equivalent form ^^'^^'^j ^ m'(-^' ^' -^i' -^2) using either hyperspherical harmonics ^2j,m,m'('^) /^) 7) 
or (complex conjugate) rotation matrix elements -D^_^/(a,/3,7) expressed in terms of Euler 
angles a, /3, 7 [85] . 

1/2 



Here j runs over 0, ^, 1, |, • • •, while m and ?n' take values in the set — 1, . . . , — j}- On 
setting 



i/3, Ai = i(7-a), A2 = 5(7 + a), pi = m - m' , p2 = -m-m' 



or, conversely, 

a = A2-Ai, /3 = 26*, 7 = A2 + Ai, m = ^{pi-p2), m = -^{pi + P2), 

and on using the relation between rotation functions d;^ and Jacobi polynomials [85], we 

indeed get 

*:;:W„p.(^,^,Ai,A2) = (-l)5(l^'^l+^^)+b^l^r°-.^^„,(fl,^,Ai,A2) (4.1) 

with 

1 /2 

^n:^,,™,^'(^,^,A„A2) = (-1)^-' ( ,2{'' + 2T+l)! ) "^-^'^ 

-m' (20)e-'('"+'"')^i gi(m-m')A2 ^ 

/:(f;)(z) = z^■42/+l)(z)e-i^ z = R\ (4.2) 

The advantage of this new form is that the wavefunctions ^nr%,pi,P2('^' ^' ^ii ^2), which were 
classified according to 

su(4) D so(4) 
[N] (2j) 

with = 2nr + 2j and 2j = 2n + \pi\ + \p2\, now turn out to be explicitly reduced with respect 
to 

su(4) D so(4) ~ su(2) © su(2) D u(l) © u(l) , . 

[N] (2j) =^ [i] © [j] H©[mr 

This will allow us to use the full machinery of angular momentum theory for determining the 
explicit action of the symmetry and dynamical algebra generators on wavefunctions. 

The two su(2) algebras appearing in chain ()4.3p are generated by Jj and Ki, i = 1,2,3, 
defined in terms of L^^, = 1, 2, 3, 4, (see equation (|2.16p ) by 

Ji = \ {\^ijkLjk — Li4) , Ki = ^ (^^eijkLjf: + Li4j , (4.4) 

where i, j, k run over 1, 2, 3 and eijk is the antisymmetric tensor. The operators Ji and Ki 
satisfy the relations 

[Ji, Jj] = icijkJk, [Ki, Kj] = icijkKk, [Ji, Kj] = 0, 
4 = Ji, Kj = K,. 
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Instead of the Cartesian components of J and K, we may use alternatively Jq = J3, J± = Ji=biJ2, 
Kq = i^3, K± = Ki lb 1^2, with Jq and Kq generating the two u(l) subalgebras in (|4.3|) . 

The differential operator form of Jq, J±, Kq, and K± can be obtained by combining equations 
(i2TD . (fZT3D . (EH]), (|2J6]1 . (I2J9]1 . and g3D and is given by 

Jo = ^{dx,-dx,), J± = ie^'(^i^^2)[±a9-i(cot0aA, +tan05Aj], 
i^o = '^{dx,+dx,), K± = ie^'(^i+^2)[T5e + i(cot0aAi -tan^aAj]. 

From some differential equation relations satisfied by rotation functions d-'^_^,{26) [86], it is 
then easy to check that 

'^O^nr,j,m,m' ~ Ur ,j,m,m' i nr,j,m,m' ~ [0 "F i m + 1)] ^ ^ rir ,j ,m±l ,m' i 

T/- vjvosc m'V[f°^'^ K V[f°^'^ \( i =n m'\(n -I- rn' -I- 1 "ll -'^/^ 11/°^*^ 

-"■0 *7ir j',m,m' ~ '''' ^ nr,j,m,m' ^ nr,j,m,m' ~ LU "T '"' jU ^ '"' ^ nr,j,m,m'±l^ 

which proves the above-mentioned result. 

(s t) 

It is now convenient to rewrite all operators of physical interest as components T^^t , cr = s, 
s — 1,..., —s, T = t,t — l, . . . , —t, of irreducible tensors of rank (s, t) with respect to su(2)©su(2). 
These must satisfy commutation relations of the type 

[Jo,T(f] = ^T(f, [J±,nf] = [{s T a){s ± a + l)]^/^rig,, 
[iro,r(f;*)] =rr(f/), [Ki,T!,f] = [{1^ r)it ± r + Ip^rj^fl,. 

Since the bosonic creation and annihilation operators serve as building blocks for the con- 
struction of other operators, let us start with them. The creation operators can be written as 
components AI-,t, cr, t = \, of an irreducible tensor of rank (|, 

The same is true for the annihilation operators, the corresponding components being given by 

= (-l)^^^-^(^Vr)^ = (4-6) 

On coupling an operator A'^ with an operator A according to 

[A^ X AY;'^ = Y.(\<y\\'y - <y'\sa){\T' ,\t - T'\tT)AlyA^^,.^r~r', 
cr',r' 

where (, |) denotes an su(2) Wigner coefficient [55], we obtain the su(4) symmetry algebra 
generators classified with respect to chain (14. 3p . These include 



J„=[A^xAy;l, = [X^ X AY^^^, c7,r = +1,0,-1, 

with J-i-i = =pJ-|-/-v/2 and K±i = ^K±l as well as the nine components of an irreducible 
tensor of rank (1, 1), 

r.,. = [^^x^]^;^, a,T = +1,0,-1. (4.7) 
The latter may be written as 

r±i,±i = -i(Tn ± 2iTi2 - T22), r±i,o = ^(±^13 - iTi4 + iT23 ± T^^), 

7±l,=Fl = -3(733 + 2ir34 - r44), 7o,±l = ^(='=^13 + i^l4 + i723 + ^24), 
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70,0 — I (^11 + T22) — — 1(?33 + T44) 

in terms of the operators T^jy, defined in ()2.17p . Observe that the u(4) first-order Casimir 
operator (I2.15P is, up to some constants, the scalar that can be obtained in such a couphng 
procedure, 

[^tx^];;; = i(^:-2). 

Similarly, the coupling of two operators provides us with the raising operators belonging 
to sp(8,]R), 

pt = . _4t = _2[^t X aX''o^ Vl, = [A^ X ^t]M, ^^r = +1,0, -1, 
or, in detail, 

= dI + dI, + dI + dI 

and 



Ko = k-Dli-Dl2 + Dl + Dl) 



2 

in terms of dJ^i, defined in ()2.18p . The corresponding lowering operators are then 
V = (Pt)t, p^^^ = (-l)-+-(Pl,_J^ a,T = +1,0,-1. 

It is now straightforward to determine the action of AI-t on the wavefunctions ^ ^/(R, 0, 
Ai, A2). Application of the Wigner-Eckart theorem with respect to su(2)©su(2) [85] indeed leads 
to the relation 

■^i,rKT,j,m,m' = Yl H W^^l^^r j) {j m, i a\j' m + a) {j m', i r|/ m' + r> 

^ ^r?^,j',m+cT,m'+T^ (4-8) 

where (n^, j'||^^||nr, j) denotes a reduced matrix element, the summation over j' runs over 
j + \i j — \^ and n'r is determined by the selection rule n'^ + j' = nr + j + \ implying that 
n'j. = HrjUr + 1, respectively. To calculate the two independent reduced matrix elements, it is 
enough to consider equation (]4.8p for the special case m = m! = j and to use the differential 



operator form of A^ i _i_ i ; 



A^ - -e^'^i 
-^±i±i " 2^ 



i ( sin 6dp> + — cos Bdg] ± — esc 6d\-^ — '\R sin ( 
KIR, 



following from ()2.ip . (I2.13|) . ()2.19p . and (14. 5|) . Simple properties of the rotation function 
<^rn -m' (2^) ^^"^ ^he Lagucrre polynomial L^^'"''^^ {z) then lead to the results 

■(2i + l)K + 2j+2)Ni/2 



{Ur.j + i||^'^||7V,i) = i 



2j + 2 
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(n.,+ 1,; - l||^t||„.„j> = -i ( Pi + l)K + l) y^\ „ 

The operators Aa,T satisfy an equation similar to (j4.8p with (n'^, j'||^"^||nr, j) replaced by 
(n^, j'||^||nr., j) and n'^ = rir — l,nr for f = j + — ^, respectively. The corresponding 
reduced matrix elements can be directly calculated from the relation 

(^r;/||-^lhr, j) = '^T^^^{nr,j\\A^\\n^,j)*, (4.10) 

which is a direct consequence of (j4.6p . 

For the su(4) generators that do not belong to so(4), we get the equation 

Ta,r^nr,j^^,m' = ^{'T''r^j'fhr,j){j m, 1 a\j' m + a){j m' , 1 t| / m' + r) 

^ ^li'^,j',m+cT,m'+T-: (4-11) 

where j' = j + 1, j, j — 1 and = — 1, n^, n,, + 1, respectively. Equation (|4.7p and the coupling 
law for reduced matrix elements |85j enable us to determine 

1/2 



/ T •^ill-^ll A f {2j + l)nr{nr + 2j + 2) \ 
{ur - l,j + l||r||n„, j> = - ^ J 

{nr,j\\T\\nr,j) = Ur+j + 1, 

(n, + u - i||r||,v,i> = - (p™i±iK!v±^) 



1/2 



from (gj]) and (liJOl) . 

The operators 'Dt,T and Po-jT satisfy a relation similar to (14. lip with 



(n„i + l||l,t||„.,i> = - ( fe- + l)K + 2,- + 2)(„. + 2,+3) ^ 



1/2 



and (n^, / 1 I'D 1 1 n^jj) obtained from these as in (j4.10p . 
Finally, with the equations 

T^Ki:j,m,m' = -2K(n,, + 2j + l)]'/'KT-l,j,n.,m', 

the action of the harmonic oscillator symmetry and dynamical algebra generators on ^^^j ^ 
9,Xi,X2) is completely determined. 

4.2 SW system potential and dynamical potential algebras 

In two dimensions, equations (12. Sp and (13. Sp simply lead to 
/c2 = 6(6-1), A;2 = a(a-1). 
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In the following, it will prove convenient to use a and b instead of ki and k2- Up to the same 
phase factor as that occurring in (|4.ip . the extended SW Hamiltonian wavefunctions (|3.9|) can 
then be rewritten a^ 



= cos'^0sin''(/)Pi""^'''"^Vcos2,^), 



2„+a+6+i^^, „i (2n + a + + a + 6 - 1)! ' 



1/2 



(n^ + 2n + a + 6)! (n + a - i) ! (n + 6 - i) ! I 
where n, a, 6 are related to j, m, m' used in (j4.2p through the relation^ 

j = n + i(a + 6 — 1), m = i(a — 6), m' = — ^(a + 6 — 1), 
or, conversely, 



(4.12) 



a = m — m + ^, 



b = —m — m' + ^, n = j + n^' ■ 



The generators of the su(4) potential algebra, as well as those of the w(4)©ssp(8, M) dynamical 
potential algebra, can be directly obtained by performing transformation p.4p on the operators 
of Section 14.11 We get for instanc^ 



J+ = -e 
^ 2 



=Fi(Ai-A2) 



=Fi(Ai+A2) 



K± = -e 



T — ^ o-2iAi 

- 4^" 



±d^ - cot i(9Ai ± o - tan (/) ic^Aa T 



=F(90 + cot |^i(9Ai ^ 9 j ~ ^ y^^-^ ^ 
— sin^ (/)(9,^ sin </> cos </>i9,^0 H — -5^Ai 2 "^^^^ "^"^l 



H — cot (j)d'^Xi "I — 2 ^^^^ 'i^^Ai "I — {\ + s\v? (j))dr -\ — 2 (cot + sin cos 



^2 CSC </'5ai - ^ 



csc^ (t!) + tj^r^ sin^ 1 



V 



A' 

^2 '^2 



i ( sin (pdr H — cos 09,*) ± - esc (^(9ai — esc — ir sin ( 

\ r / r 2r 



1 



-2iAi 



sm 



,2 i o2 



sm (/) cos (pO~^ 



,2 , 2iq2 



1 2 o2 o2 

COS (pa J, H — ^ cot 



1 



+ _csc^</.5^^ + 



— (1 + sin (/)) + 2ti;r sin 
r 



dr + 2 



(cot + sin (f) cos ( 



+ 00 sin ffi) cos I 



— i ( — csc^ + 2a; j 9ai — c^c^ — w^r^ sin^ — w ^ , 



*It is worth observing here that integer or half-integer values of j, m, and m' are related to integer values of n 
and half-integer ones of a and b. The results for matrix elements of potential and dynamical potential algebra 
generators are only valid for such a and b (see footnote [3]) , although those for wavefunctions are not restricted to 
these values provided factorials are replaced by gamma functions. 

^Equation (|4.12p is valid for positive pi and p2, corresponding to physical wavefunctions (see discussion at the 
end of Section [3]). 

®For simplicity's sake, we denote both types of operators by the same symbols. 
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= J- (-H - 2ujrdr + 2a; - 2a;) . 
2a; ^ 

It is also straightforward to derive their matrix elements from the results of Section 14.11 and 
equation (|4.12|) . We list them below: 

Jo'^nr,n,a,b = ^(o - b)^nr,n,a,b, Ko^nr,n,a,b = -^{a + b - l)^nr,n,a,b, 

J+^nr,n,a,b = [{u + a + \) [n + h - \)] ^'„,,„,a+l,b-l, 

J-^nr,n,a,b = [{u + a - \) [n + h + \)] ^'n„n,a-l,6+l> 

K+^n^^n,a,b =[{n + l){n + a + h- l)\^/'^^nr ,n+l,a— 1 ) 

K_^rir,n,a,b = [n{n + a + 6)]^/^^'„^,„_i,a+l,6+l, 
n+r+1 

Ta,T^nr,n,a,b= ^ t„/ (n^ , 2n + a + 6) 
n/=n+T— 1 

X (n + i(a + 6 - 1) i(a - 6), 1 cr|n' - r + i(a + 6 - 1) i(a - 6) + a) 

X (n + i(a + 6 - 1) - i(a + 6 - 1), 1 T\n - r + i(a + 6 - 1) - i(a + 6 - 1) + r) 

^ nr — {n' —n—T),n' ,a+(T—T,b—a—T ^ 
n+r+i 

•^l,T^nr,n,a,b = ^ ^n' (n^ , 2n + a + 5) 
n'=?i+T— i 

X (n + i(a + 6 - 1) i(a - b), ^ a\n' - t + ^{a + b - 1) ^{a - b) + a) 

X (n + i(a + 6 - 1) - i(a + 6 - 1), i T|n - r + i(a + 6 - 1) - ^(a + 6 - 1) + r) 



X ^ 



rir — (n'— n— r)+i n',a+(T— T,f)— cr— T' 



2 ' 

n+r+1 



^i,T*n,,n,a,6 = ^ ^n' (n^, 2n + a + b) 
n'=n+T—l 

X (n + i(a + 6 - 1) i(a - 5), 1 cr|n' - r + i(a + 6 - 1) i(a - 6) + cr> 

X (n + i(a + 6 - 1) - i(a + 6 - 1), 1 T\n - r + i(a + 6 - 1) - i(a + 6 - 1) + r) 



X vl/ 



rir — (n'— n— r) + l,n',a+(T— T,fe— cr— r 5 



-2[K + 1)K + 2n + a + 6 + l)]^/2^'n.+l,n,a,6- 



Here 



t^' (n^, 2n + a + 6) = < 



and 



dn'{nr, 2n + o + 6) 



/ (2n+a+b)n, (n,+2n+a+b+l) \ If _ „ , ^ i 1 

\^ 2n+a+6+2 ^ U n — -|- T i- i, 

?ir + ^^ + ^(o + & + 1) if n' = n + r, 

{2n+a+b){nr+l){nr+2n+a+b) \ "^/^ 



2?i+a+6-2 



if n' = n + r — 1 , 



o„' (n,. , 2n + a + 6) = < 



. f {2n+a+b){nr+2n+a+b+l) \ ^1'^ 



2n+a+b+l 



_. / (2n+a+b)K + l) y/^ 
\^ 2n+a+fe-l y 



j ifn' = n + T + ^, 
if n' = n + T — ^ , 



(2n+a+6)(nr+2n+a+fe+l)(nr+2n+a+6+2)\ -"-/^ -f / _ , , i 



2n+a+fe+2 

[(n,. + l)(n^ + 2n + a + 6 + 1)]^/^ 

(2n+a+b)(n^+l)(nr+2) \ 



2n+a+6-2 



if n' = n + T, 
if n' = n + T — 1. 
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From these results, we conclude that the potential algebra generators produce transitions 
between levels belonging to spectra of Hamiltonians characterized by parameters (a, b), (a it 1, 
6^1), (a ± 1, 6 lb 1), (a ± 2, b), and (a, b it 2). For the dynamical potential algebra generators, 
the same Hamiltonians are involved together with those associated with (a it 1, 6) and (a, 6 it 1). 



5 Conclusion 



In the present paper, we have re-examined the L'-dimensional SW system, which may be con- 
sidered as the archetype of D-dimensional superintegrable system. We have completed Evans 
previous algebraic study, wherein its symmetry and dynamical algebras had been determined, 
by constructing its potential and dynamical potential algebras. 

In our approach based on the use of hyperspherical coordinates in the L'-dimensional space 
and on the introduction of D auxiliary continuous variables, the SW system has been obtained 
by reducing a 2Z)-dimensional harmonic oscillator Hamiltonian. The su(2D) symmetry and 
w(2L')©sSp(4Z),M) dynamical algebras of the latter have then been transformed into correspon- 
ding potential and dynamical potential algebras for the former. Finally, the two-dimensional 
case has been studied in the fullest detail. 

Possible connections with other approaches currently used in connection with the SW system 
or, more generally, superintegrable systems, such as supersymmetry [521E3], path integrals [87] . 
coherent states [88], and deformations \89\ 190]. might be interesting topics for future investiga- 
tion. 

A Wavefunctions of the 2D-dimensional harmonic oscillator 

The purpose of this appendix is to derive the explicit form of the harmonic oscillator wavefunc- 
tions (fZ^ . 

On inserting ()2.3p in the Schrodinger equation ()2.2p . the latter separates into D — 1 angular 
equations 

\-dl - [{2D -2v-l) cot 9, - i3ii9,]de,^ + + - cA @,{9,) = 0, 

[ sm^ 9y cos^ 9u J 

u = l,2,...,D-l, (A.l) 
and a radial equation 

-4 - + § + - ^°") ^(^) = 0- 

Here Ci, C2, . . . , Cz)_i are D — 1 separation constants, while C^) is defined by 

Cd=pI. (A.3) 

In the following, we are going to show that there does exist a solution to the whole set of D 
equations ()A.ip and ()A.2p such that all the separation constants C,^, v = 1,2,...,!) — 1, are 
nonnegative. 

Let us start by solving the angular equation ()A.1|) corresponding to the variable 9y in terms 
of and Cyj^\. The ansatz 

Q^{9^) = {cos9^f''-^sm9^)''''-^F^{u^), = cos^ 9^, 

transforms it into the hypergeometric differential equation [84] 

{u^{l - Uu)dl^ + [7 - (a + ^ + l)u^]du, - al3] F^{u^) = (A.4) 
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provided we choose the constants and 6^ in such a way that 

{a,-^f=pl, {K-^){b, + 2D-2i^-^) =C,+i. (A.5) 
In (|A.4p . a, 13, and 7 are given by 

a = ^{a^ + bu + D-iy-l + A,,), p = \{ay + K + D - v - I - /^^), 7 = 0^ + !, 

where 

A, = ^{D - vY + C,. (A.6) 

There are altogether four solutions to the two quadratic equations (lA.Sh . which may be written 
as 

a^ = \ + e\py\, K = -{D-v-l)+eAu+i, e,e' = ±l. (A.7) 

Consequently, we get 

a = i(l + e\p^\ + e'A^+i + A^), P = \{l + e\py\ + e'A^+i - A^), 

7 = l + e|pj.|. (A.8) 

The general solution of the differential equation ()A.4p may be written down as 

F^K) = A 2Fi(a, /?; 7; Uy) + Bul^^ 2^^1(0 - 7 + 1, /? - 7 + 1; 2 - 7; ""i^), 

where A and i? are two constants to be determined so that the angular function Qu{Ou) be 
physically acceptable, i.e., vanish for — )■ and 9y ^ ^. On considering the four possibilities 
for the pair (e, e') in (IA.7|) and (lA.Sp successively, we arrive at a single solution corresponding 
either to e = +1, e' = +1, S = 0, /3 = -n^ {n^ e N) or to e = -1, e' = +1, ^ = 0, 
^ — 7 + 1 = — rijy (rijy G N). It can be expressed in terms of a Jacobi polynomial [84] as in 
equation ()2.7p . where 

au = \Pu\ + h b,y = -{D - u - + A^+i, = 0,1,2,..., (A.9) 

while the separation constant must satisfy the equation 

= 2n^ + \p^\ + A^+i + 1 (A.IO) 

with Ai, defined in ()A.6p . 

To obtain a solution to the whole set oi D — 1 angular equations (jA.ip . as expressed in 
equation (j2.6p . it only remains to solve the recursion relation (jA.lOp for A^, with the starting 
value Ad = \pd\ corresponding to ()A.3p . The results for Ajy and Cu, u = 1,2, . . . , D — 1, read 

Ai, = 2ny + 2ny+i H + 2nD-i + \pv\ + \py+i\ H h \pd\ + D - ly 

and 

Cu = {2ny + 2n^+i H h 2nD^i + \pu\ + \Pu+i\ H h \pd\) 

X (2n,. + 2n^+i + ■■■ + 2nD-i + \Pu\ + \Pu+i\ H + \pd\ + 2D- 2v), (A. 11) 

respectively. As a consequence, 61, in (|A.9p can be rewritten as in equation (|2.8p . This completes 
the proof of equations ()2.6p to ()2.8p . 

Turning now ourselves to the radial equation (lA.2p . we note from (lA.lip that Ci can be 
written as 

Ci = Aj{j + D-1) 

in terms of j defined in ()2.10p . Finally, it is straightforward to show that the physically accept- 
able solutions vanishing for r (or z) going to zero and infinity are given by (12. 9p and correspond 
to the eigenvalues (j2.1ip . 
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